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Abstract We provide a 2.5-dimensional solution to a complete set of viscous hydrody- 
namical equations describing accretion-induced outflow and then plausible jet around black 
holes/compact objects. We prescribe a self-consistent advective disk-outflow coupling model, 
which explicitly includes the information of vertical flux. Inter-connecting dynamics of 
inflow-outflow system essentially upholds the conservation laws. We provide a set of an- 
alytical family of solutions through the self-similar approach. The flow parameters of the 
disk-outflow system depend strongly on viscosity parameter a and cooling factor /. 

1 INTRODUCTION 

Most extragalactic radio sources are expected to form around spinning massive black holes (Meier et al. 
2001 Meier 2002). The immense amount of matter, forming an accretion disk, is being accreted either 
from the interstellar medium or from its companion star. In these systems, the relativistic outflowing matter 
should come only from the inner regions of the accretion disk unlike stellar outflows. This is particularly 
suggestive for quasars or the micro-quasars which do not have an atmosphere of their own. Fender, Belloni 
& Gallo (2004) suggested a semi-quantitative model for the jet in black hole X-ray binaries where a cor- 
relation between the radio and the X-ray emission was estimated. Vadawale et al. (12001) established the 
X-ray and radio properties of micro-quasar GRS 1915+105. Time dependent interaction between the jet 
and the inner disk (e.g. Ueda et al. 2002 ) was evident from the observations of simultaneous X-ray/IR flares 
from a black hole/relativistic system. Rawlings & Saunders ( 119911 ) found a strong correlation between the 
narrow-line and radio luminosity in FRII type radio galaxies. This implies that the production of optical 
line emission and large-scale radio emission are intrinsically linked. Therefore, the outflows or jets are ex- 
pected to correlate with the disk controlling the accretion process, precisely the accretion dynamics around 
a central star. The jets or outflows extract matter, energy and angular momentum from the disk. 

Thus it is now clear that these two apparently dissimilar objects are related each other. In principle, 
one should study the disks and outflows leading to jets in a unified manner, which cannot be dealt as 
separate flow dynamics. However, there are few models which simultaneously study the accretion-outflow 
dynamics on the same platform. Chakrabarti & Bhaskaran (1992) attempted to correlate the collimated 
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and self-similarity in the radial direction. Blandford & Begelman d 19991 ) modified the ADAF solution, 
originally proposed by Narayan & Yi ( 1994) where the accretion flow is well below the Eddington limit, by 
including an outflow/wind which carries mass, angular momentum and energy from the accretion disk. They 
later extended their work to two-dimensional adiabatic flow (Blandford & Begelman 2004). Although a new 
branch of wind solutions was discovered, that does not include the vertical fluxes in the hydrodynamical 
equations. The ADAF model (Narayan & Yi |1994|[T99 5 ) explained the under-luminous accreting sources. 
The interesting aspect of the ADAF model is that the Bernoulli's parameter at all radii (within the acceptable 
location of the validity of the self-similar approach) is positive, which leads to conceive that the outflows 
and jets might emanate from the advective disk. Later on, stability of the solution under perturbation was 
studied with the inclusion of Coriolis force by Prasanna & Mukhopadhyay (2003 ). In recent times a few 
simulations on disk-outflow coupling have been cultivated (Nishikawa et al. 2007 McKinney & Narayan 
2007 ). However, the results are strongly dependent on the initial conditions (Ustyugova et al. 1 19991 ) and it is 
difficult to simultaneously simulate the disk and the outflow regions because the time scales of the accretion 
and outflow are in general very different. Moreover, in these simulations how the matter gets deflected from 
the equatorial plane has been studied largely in the Keplerian regime. 

In recent years, there have been a discovery of unusual class of compact sources, the ultra-luminous X- 
ray sources (ULX), in the nearby star forming galaxies (Katz l 19871 Fabbiano et al. 1989 Kaaret et al. 1200 II 
Colbert & Ptak 2002, Miller et al. 2003 , Begelman et al. 2006 ). These are optically thick, radiation pressure 
dominated systems with strong advection and the matter is strongly ejected out from the disk in the form of 
outflows/jets by strong radiation pressure. Using a slim disk model, Abramowicz et al. ( 119881 ) discovered 
a new branch of solution at a super-critical rate which is stable and optically thick. A model for super- 
critical accretion with advection was attempted by Lipunova ( 1999 ). Ohsuga et al. (2005 ) have emphasized 
the importance of advective flows in the super-Eddington, radiation pressure dominated disk with photon 
trapping. Hence, these two opposite paradigms of black hole activities reveal a profound inter-connection 
between the inflow parameters and the outflows leading to jets, especially in the advective regime, which 
the standard optically thick Keplerian disk theory fails to explain. 

In the present work, without assuming a geometrically thin disk structure, we prescribe a new model for 
the accretion-induced outflow leading to jet. We construct the inflow-outflow correlation model in a more 
self-consistent manner. The contribution of magnetic field is neglected at the first instant. The magnetic field 
is more important to explain the collimation and acceleration of jet (apart from ultra luminous sources]^ 
The present model can, not only extend our model from quasars to micro-quasars, but also to neutron 
star X-ray binaries and in general to many sources with outflows from the disk. However, to describe the 
flow dynamics and consequently outflows in protostellar objects the standard Keplerian disk model itself is 
enough. Our unification scheme is based on the fact that the astrophysical outflow and jet, its underlying 
disk and its inter-related dynamics at all scales, obey same physical laws. 

We arrange our paper in the following manner. In the next section, we formulate our model equations 
for the accretion-induced outflow. In §3, we present a complete analytical, but self-similar, solution of our 
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model. Next, we study the properties of the class of solution in §4. In §5 we end with a discussion and 
summary. 

2 DISK-OUTFLOW CORRELATION AND MODEL EQUATIONS 

We assume the disk to be steady and axisymmetric. For a generalized geometrically thick advective disk, 
we consider the r<fr— , <fiz— and rz— components of the shearing stress. The remaining stresses are believed 
to be negligible which do not significantly contribute to control the disk-outflow dynamics. The flow pa- 
rameters v r , A, v z , c s , p and P are considered to be functions of both radial and vertical coordinates, which 
are radial velocity, specific angular momentum, vertical velocity, adiabatic sound speed, mass density and 
pressure respectively. Here, throughout our calculations, we express radial and vertical coordinate in the 
unit of 2GM / c 2 , where M is mass of the central star, G is the gravitational constant and c is speed of light. 
We also express velocities in the unit of speed of light and specific angular momentum in 2GM/c. The mass 
of the disk is assumed to be much less than that of the central object, hence the disk is not self-gravitating. 
Therefore, the general disk-outflow coupled equations are given below, 
(a) Mass transfer: 
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where W rz is the rz th component of the stress tensor, when we consider the shear stress tensor is sym- 
metric (Landau & Lifshitz 1989 ) and For is radial component of the gravitational force. To understand the 
importance of the term dW rz /dz in the above equation, we compare it with dP /dr as 
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where we use a generic order of magnitude relation dA/dxj 0(A/Xj); A denotes any independent 

quantity as a function of an arbitrary coordinate variable Xj, h(r) is the disk half-thickness. Note that we do 

not identify h here as a hydrostatic scale height, instead the photospheric height where the disk is coupled 

to the corona, vt is the turbulent kinematic viscosity. With c 2 ~ Pj ' p and from eqn. ([TJ we obtain 

I v z | h 
' v r r 

However, we can write from eqn. (0 (as described below) 
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Using eqns. (0]i and ©, and assuming an isotropic distribution of turbulencq 2 ] such that v t ~ ac s h, where 
a (a < 1) is the Shakura-Sunyaev viscosity parameter (Shakura & Sunyaev 1973), eqn. (O reduces to 
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For a reasonable value of h ~ r/2 and a, the second term on the right hand side of eqn. © can be neglected. 

dP 



Thus we retain with (dW rz /dz) / 4£ ~ a 2 ; dW rz /dz can not be neglected. In order to determine W rz , 



we derive a simplified relation of W rz with W r< £, which is the r<fi th component of the stress tensor, from 
the order of magnitude analysis and obtain W rz ~ aW r< /, h/r 

As the disk has a significant radial flow, we include ram pressure along with gas pressure 
(Mukhopadhyay & Ghosh 2003) in the equations and write W r< j, = —a(P + pv 2 .). The radial momen- 
tum equation of the disk-induced outflow/jet thus reduces to 

dv r dv r A 2 _, 1 dP a 2 d ,z . 2x1 

v r ^~ + v z - 3 + F Gr + -— + — — {-(P + pv 2 )} = 0, (7) 

or oz r° p or p oz r 

where we have used the fact that for a thick disk, in general, h ~ z. 
(c) Azimuthal momentum balance: 
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The first term on the right hand side signifies the outward transport of angular momentum in the radial 
direction and the second term in the vertical direction due to the turbulent stress. 

If W r< j) dominants the angular momentum transport, with the use of mass conservation eqn. (Q~|), we 
obtain 
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If, on the other hand, W^ z dominates the angular momentum transport, then we obtain 

, , \Wa % \ 



Now comparing eqns. (O and (fTOb and with the use of mass conservation eqn. (Q]i, we can write 
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Therefore, the azimuthal equation reduces to 
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(d) Vertical momentum balance: 
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where F Gz is the vertical component of the gravitational force. As before we estimate 
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eqn. $13[ thus reduces to 



a (\ + \). For h ~ r/2 and reasonable a, the quantity is negligible. The 
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In absence of first term, the equation leads to a mean vertical outflow from the disk. On the other hand, 
if there is no outflow and jet: v z = 0, then eqn. ( TT4l > reduces to the well known hydrostatic equilibrium 
condition in the disk, from where one can calculate the hydrostatic disk-scale height, 
(e) Energy conservation: 

For an accretion-induced outflow the energy budget can be computed by 

~(r^ P ) + ^=0, (15) 
r or oz 

where T r and T z are the radial and vertical components of the total energy flux Ti given by 

(v 2 7 P \ 
Fi = pVi[— + +<Pa)- VjWij + F h (16) 
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where we neglect the molecular heat conduction as this is insignificant in accretion disks and the nuclear 
heat generation/absorption (Mukhopadhyay & Chakrabarti 2000) for mathematical simplicity. Here v 2 = 
v 2 + X 2 /r 2 + v 2 , 4>g is the gravitational potential, Fi the radiative flux from the disk surface, Wij the 
generalized stress tensor, 7 the gas constant: 4/3 < 7 < 5/3. 

Using eqns. ([1]), ©, ||S), ( fT3] > & dT6b , we obtain the disk-energy equation from ( fT3T > as 



pv ■ Vs = 



r, - 1 



dP Pdp 
dr 1 p dr 



dP Pdp 

dz 1 p dz 



= Q + -Q- = fQ+, (17) 



r 3 -i 

where s is entropy density. Here we assume the energy released Q~ due to radiative loss from the disk 
is proportional to the viscous heat generated, Q + , where / is the cooling factor incorporating any kind of 
outflow and jet which is close to and 1 for the flow with efficient and inefficient cooling respectively. 
The first and the second terms on the left hand side of the above equation represent the radial and vertical 
advection of the flow respectively, where we define (Cox & Giuli |1968l Mukhopadhyay & Ghosh |2003T > 
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with (3 = |^5§, the ratio of gas pressure to total pressure which is close to for extreme radiation dominated 
flow (7 = 4/3) and to 1 for extreme gas dominated flow (7 = 5/3), a is Stefan constant, m p is mass of the 
proton, T is proton temperature, ks is the Boltzmann constant, p is average molecular weight. 
In eqn. ( [T7| >, Q + = W^/rft, f]t is the coefficient of turbulent viscosity. Thus for our case 

Q + = ^{W^ + W 2 Z + W 2 Z ). (19) 

As before, it can easily be shown that the contribution of W rz is much less than that due to W r( p 
(Wrz/W 2 ^ ~ a 2 [^n + ^r] ^ . W<p z contributes to the additional viscous heating in a geometrically thick 
advective disk with vertical outflow. Using mixed shear stress formalism (Chakrabarti 1996) and approxi- 
mating Wfa in terms of W ri p as given by eqn. (fill , eqn. ( fTTl ) reduces to 

Vr\dP r Pdp] v z \dP r Pdp] 2 lfdX A M zdX\ 
Ti — — + Ti — — = -fa (P + pv r - —-2- + -— .(20) 
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3 SOLUTION AND SELF-SIMILARITY 



We follow the self-similar approach to solve the equations in obtaining the class of solutions. For the present 
purpose we seek for a generalized self-similar solution, unlike the previous case (Narayan & Yi 1994), 
where variation of the flow parameters as functions of vertical coordinate along with radial coordinate has 
been invoked for a coupled set of disk-outflow equations of the form 



v r (r,z) = v r0 r u z a , X(r,z) = X r v z , v z (r,z) = v z0 r w z , c s (r,z) = c s0 r 9 z s , 



(21) 



where v r o, Xq, v z q and c s q are the dimensionless coefficients which will be evaluated from the conservation 
equations. We determine the exponents u, a, v, b, w, d, g, s by self comparison of various terms in the equa- 
tions. Assuming the flow to be polytropic, as most likely the disk is, we consider the adiabatic equation of 
state as P = kp 1 , where 7 = 1 + 1/n, n is the polytropic index of the flow, while the adiabatic sound speed 



We propose a generalized gravitational potential 0g( 7 ' 5 z ) = —{r — ~i^f2 r ~ 3z2 ) zk > where the index 
k induces the variation along z-axis which we determine self-consistently. When the disk does not have 
strong outflow and jet k = 0, and <\>q reduces to conventional Newtonian potential upto the second order in 

(z/r). 

Substituting the solutions from eqn. (l2TT i in eqns. (fTJ and (0 and comparing the exponents of r and z 
we obtain w = u — l,o = s, d = a + 1, u = —1/2, v = 1/2, g = —1/2, b = a, k = 2a. 
Equations ([TJ, Q, ( fl2l l and ( f20b can now be written, with the use of eqn. (1211 ). respectively 
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Solving eqns. d22T> - (T23T> we compute the coefficients of eqn. (f2TJ 
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1 
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where Q = X + 7 V 2 , X is given by 
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and A, B, V, £, H and K, are represented as 
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The generalized Bernoulli equation is then 
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where S^; is the Bernoulli constant. 

The above solutions can explain both super-critical and sub-critical accretion flows, where the flow is 
more likely to be strongly advective with strong possibility of the outflow and jet. Super-critical accretion, of 
the order of M ^ (1CP 3 — 1CP 6 )M© /yr, corresponds to high luminosity sources with mass of the central 
star M ~ IOMq. In this case, the flow is expected to be radiation pressure dominated with maximum 
physically plausible 7 is 1.444, corresponding to P r ~ P g . 

To determine the exponent a, we vertically integrate eqn. ( TBi i from —h to +h after substituting the 
solutions given by eqns. d2TI ). (t26b-(|29b. As the outflow is not likely to emanate from the equatorial plane, 
the solution looses its relevance there because the torque due to W^, z exerted on the matter is zero. However, 
from a certain finite height ho, they are relevant describing a disk-outflow system. We consider h/r ~ t < 1 
and ho/r ~ to <C 1, where t and to are kept constant throughout our analysis. The realistic flow, when v r < 
and v z > 0, demands that a cannot be positive. This helps us to fix the boundary condition of the outflow in 
the vertical direction. We demand a situation for which to is least to yield a physically realistic a. For super- 
critical flows exhibiting high luminosity sources, using eqn. ( T22l i we obtain a most physically acceptable 
solution for a given by a ~ —(2/11 + e) for an appropriate to ~ 0.02 corresponding to a reasonable 
t ~ 0.5, when e is a very small number ^ 10~ 3 . For sub-critical flows exhibiting under-luminous sources, 
on the other hand, the highly sub-critical mass accretion rate M < (10~ 10 — l0~ 12 )MQ/yr or M < 
(10~ D — l0~ 7 )M Q /yr corresponding to black holes of mass M ~ 10M Q orM ~ 10 6 M Q respectively, 
for which P g 3> P r and 7 5/3. With a similar argument as above we obtain here a ~ — (1/4 + e) for an 
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4 PROPERTIES OF SELF-SIMILAR SOLUTIONS 



Here we study the properties of our solution for cases of under-luminous and highly luminous sources 
at sub-critical and super-critical accretion rates respectively. The standard model of Shakura & Sunyaev 
( 1 19731 1 is ineffective to describe these two cases of the geometrically thick advective disks having a sub- 
stantial outflow. We describe a typical set of solutions for the accretion-induced outflow using the flow 
parameters obtained in the last section in these two opposite paradigms. A detailed family of solutions and 



their observational implications will be discussed elsewhere (Ghosh et al. in prep. 

4.1 Super-critical accretion regime 

Let us consider a case where radiation pressure P r dominates over gas pressure P g for the flow with high 
Eddington-accretion rate. We choose j ~ 1.4 corresponding to j3 ~ 1/3, appropriate for the above class 
of flow. The flow is radiation trapped, optically thick and hot. Figure 1 describes the variations of flow 
parameters as functions of radial and vertical coordinates for two values of / at a typical a. We see that v r , 
v z , A and c s fall off rapidly with the increase of / at a fixed z. With the decrease of /, the outflow velocity 
v z increases rapidly due to strong radiation pressure which blows up the matter as shown in Fig. lc. The 
disk gets possibly truncated due to strong outflow having both radiation and gas at a region around r ~ 17 
for / - 0.4, and r ~ 10 for / ~ 0.7. 

In general, an increase of / leads to the inefficient cooling. This renders the disk to be puffed-up and 
more quasi-spherical. As a result, the disk angular momentum decreases due to its extraction by the out- 
flow/jet (see Figs. lb,c). However, at higher / (> 0.5), the system becomes radiatively very inefficient, 
which may result in the decrease of the possible outflux with an increase of / rendering an increase of the 
flow angular momentum. At low a (~ 0.01) when the residence time of the infalling matter in the disk 
is high, angular momentum of the system is such that the disk becomes centrifugally dominated. At this 
stage, with the decrease of /, angular momentum may increase resulting in the radial and vertical velocity 
of the flow to enhance significantly in order to overcome the strong centrifugal barrier. The outflow is then 
centrifugally dominated. 

The Bernoulli's number in Fig. 3 is similar to that of the velocity profiles in Fig. 1. Be is always positive 
and high at low /, which indicates the plausibility of the outflow to be very strong, and falls off rapidly with 
r. The probability of the outflow and jet is low at high /. With an increase of z, Be initially decreases. This 
is due to the fact that the first term of the potential <j)Q, which dominates at small z, is attractive in nature. 
Then Be gets a kick as the repulsive part of 4>q dominates with the increase in z. 

4.2 Sub-critical accretion regime 

For highly sub-critical accretion flows, which are associated with very low density plasma, the possibility 
of transfer of viscous energy from ions to electrons due to the Coulomb collisions is very negligible. This 
results in a gas pressure dominated geometrically thick accretion disk. The flows have strong advection due 
to inefficient cooling and are optically thin (Narayan & Yi 119941 [1995 ). To analyse our result we choose 
7 = 1.6 which corresponds to j3 ~ 0.89 and / = 0.9. Figure 2 shows the profiles of the flow parameters 
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Fig. 1 Variation of (a) velocity, (b) specific angular momentum, (c) vertical velocity, (d) sound 
speed, as functions of radial and vertical coordinates for super-Eddington accretion flows. Solid 
and dashed sheets are for / = 0.4, 0.7 respectively. Other parameters are a — 0.05, 7 ~ 1.4 and 
corresponding (i ~ 0.3. 

signify that the magnitudes of v r and v z are much less compared to that in the super-critical accretion flows. 
In low mass accretion flows, the disk may get truncated at much nearer to the central star. 

5 SUMMARY 

We have presented a self-consistent model of accretion-induced outflow and then jet. We have established 
our model equations in a more general way, than done earlier, without making any hypothesis, and without 
restricting ourselves to the Keplerian geometry. Our equations uphold the conservation laws as the outflows 
and jets extract matter, energy and angular momentum from the inflowing matter. In its analytical self- 
similar form, it is more easy to analyse and study the family of solutions (with variation of a&/) and to 
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Fig. 2 Same as in Fig. 1, but for sub-Eddington accretion flows. Other parameters are a — 0.05, 
/ = 0.9, 7 ~ 1.6 and corresponding j3 ~ 0.89. 

we have kept is to ignore the importance of magnetic field in the disk-outflow system. While not including 
magnetic field is an assumption, the outflows and then jets in ULX are expected to emerge due to strong 
radiation pressure. Therefore, the collimation of jet in ULX might not be magnetically linked (Jaroszynski 
& Abramowicz 119801 Fabrika 2004). Therefore, for ULX and highly luminous AGN, the assumption of 
neglecting magnetic effects could be quite appropriate. We also do not aspire to describe the mechanism 
for formation of jets, for that the inclusion of magnetic terms might be mandatory, but try to understand the 
accretion flow dynamics with the inclusion of the vertical flow. Moreover, to include magnetic field, solve 
the equations, and obtain the solution in its present form is beyond the scope. 

The new insights that we have provided in the work are: 
1) We have studied the complete set of axisymmetric Navier-Stokes equations for accretion-induced out- 
flow analytically. 
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Fig. 3 (a) Variation of Bernoulli's constant for high mass accretion flows. Solid and dashed 
sheets are for / = 0.4,0.7 respectively. Other parameters are a = 0.05, 7 ~ 1.4 and corre- 
sponding ft ~ 0.3. (b) Same as in (a), but for low mass accretion flows, a = 0.05, / = 0.9, 
7 = 1.6 and corresponding (3 ~ 0.89. 

model, thus invoking a 2.5-dimensional accretion flow. 

3) We do not assume hydrostatic equilibrium. 

4) We have explicitly included <pz— and rz— components of stress tensor apart from the usual r<j>- 
component in order to include outflow dynamics into the disk. 

5) All the flow parameters are considered to be functions of both r and z coordinates. We explicitly have 
shown, by order of magnitude analysis, which terms are relevant and which others can be discarded. 

Two extreme cases of the geometrically thick advective accretion disk consisting of super-critical and 
high sub-critical accretion flows have been studied. It shows that the dynamics of the system depends 
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to that of super-critical flows. Although we have made a self-similar analytical study, it exhibits some 
reasonable features in understanding the dynamics of the accretion-induced outflow. 
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